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ABSTRACT 

Stationary, spherical accretion onto an unmagnetized neutron star is here 
reconsidered on the wake of the seminal paper by Zel'dovich & Shakura 
(1969). It is found that new "hot" solutions may exist for a wide range of 
luminosities. These solutions are characterized by a high temperature, 10 9 -hl0 11 
K, and arise from a stationary equilibrium model where the dominant radiative 
mechanisms are multiple Compton scattering and bremsstrahlung emission. For 
low luminosities, ^ 10 -2 Le, only the "cold" (a la Zel'dovich and Shakura) 
solution is present. 
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I. INTRODUCTION 

Even before the observational evidence that Galactic X-ray sources are mostly 
binary systems containing an accreting neutron star, Zel'dovich & Shakura 
(1969, ZS in the following) studied in some detail the spectrum of radiation 
produced by stationary, spherical accretion onto an unmagnetized neutron star 
and compared their results with the (poor) data available at the time for Sco 
X-l. The pioneering paper of ZS shows that the resulting spectrum depends 
on two parameters, the accretion rate (luminosity) and the penetration length 
of the accreting ions in the outermost neutron star layers. The outcome can 
be described, in essence, as a black body with a high energy tail due to the 
Compton heating of thermal photons in the hot, external part of the atmosphere 
surrounding the neutron star. 

ZS's analytical work was pushed further by Alme & Wilson (1973), making 
use of numerical methods. Shapiro & Salpeter (1973) adopted essentially 
ZS solution at the inner boundary and explored the possibility that in 
the surrounding region a shock is formed, which may modify the resulting 
spectrum. While models of spherical accretion evolved substantially since then, 
considering, for instance, the role of nuclear reactions induced in the crust by 
the bombardment of accreting ions (see e.g. Bildsten, Salpeter & Wasserman 
1992), the basic picture proposed by ZS has been maintained. 

In this letter we reconsider the issue, adopting a set of equations for the 
neutron star atmosphere which essentially coincides with those of ZS. We show 
that ZS's solution for the temperature profile is not unique: for large enough 
accretion rates, another solution, at considerably higher temperatures, may 
exist. The presence of these solutions per se, and the transition between "hot" 
and "cold" solutions can have interesting astrophysical consequences. 
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II. THE MODEL 

In order to explore the possible existence of high temperature solutions, we use 
a very simplified model which follows closely the original work of Zel'dovich 
& Shakura. In particular we assume that the accreting flow impinges onto a 
spherically-symmetric, static atmosphere which surrounds the neutron star, and 
is decelerated as it penetrates into the atmosphere. The kinetic energy released 
by the incident protons goes mainly into electron thermal energy and is then 
re-emitted as free-free radiation. The details of the flow braking fall within the 
domain of plasma physics and are still far from a thorough understanding. For 
this reason, following ZS, we treat the total column density of the atmosphere 
required to stop the incoming beam, y , as a free parameter and the column 
density 

/•OO 

y= QdR (1) 



JR 

is used as the independent coordinate in place of the radial distance R; here g 
is the matter density. 

The heat injected by the infalling protons per unit time and mass in the 
atmosphere is assumed to be constant and is related to the total luminosity 
observed at infinity by 

Lc 



W = - ^ - ■ (2) 



In the inner region where y > yo, W = 0; in all our models 5 < yo < 20 . If 
the flow velocity v exceeds the electron thermal velocity, v t h, these values of 
the proton penetration length are appropriate to describe the stopping of the 
incoming proton beam in a hydrogen plasma where only Coulomb interactions 
take place (e.g. Alme and Wilson 1973). On the contrary, if v& v t h, the 
proton stopping through repeated Coulomb scatterings is less effective and other 
collective processes (e.g. plasma oscillations) need to be considered to keep yo 
within this interval. 

The transfer of radiation in the atmosphere is governed by the equations 
for the radiative luminosity L and the radiation energy density U which, in 
spherical symmetry and using the Eddington approximation, can be written as 

^ = -4nR 2 W (3) 
dy 
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The only radiative processes taken into account are scattering and 
bremsstrahlung, so that the flux mean opacity can be conveniently expressed 
as 

«i = K es + 6.4 x 10 22 gT~ 7/2 cm 2 g" 1 

where K es = 0.4 cm 2 g -1 and T is the gas temperature. Since, as numerical 
models show, the atmosphere does not expand considerably even for high 
temperature solutions, the Eddington approximation is reasonable. The 
appropriate boundary condition for the radiation field at the outer edge of the 
non-illuminated medium is U = L 00 /2irR 2 c. The inner boundary condition, at 
Din ^> VOi is fixed by the requirement that all the observed radiative flux must 
be generated within the atmosphere, that is, L = 0. This is the same condition 
used by ZS and is appropriate if the atmosphere is effectively thick close to yi n . 

The runs of pressure, P, and temperature are obtained from the hydrostatic 
balance and radiative energy equilibrium 

dP _ GM* 
~dy ~ E? 

where M* is the mass of the neutron star. The matter density g is calculated 
from the perfect gas equation of state assuming that the material is completely 
ionized hydrogen. Since we study the general properties of spherical accretion 
onto neutron stars for luminosities below the Eddington limit, in equation (5) 
we neglected the radiative force and the ram pressure exerted by the incoming 
protons, which is typically two orders of magnitude below the gravitational 
term. In equation (6) Bp is the Planck mean opacity and T 7 is the radiation 
temperature which is defined as the mean photon energy. In general T 7 can be 
computed only solving the full frequency-dependent transfer problem and will 
depend on y. In ZS, T 7 was taken equal to \U{y)/aY^, which is appropriate in 
LTE. Being interested also in solutions in which multiple Compton scattering 
becomes important, we derive the radiation temperature from the equation (see 
Wandel, Yahil & Milgrom 1984; Park & Ostriker 1989) 



(5) 



y dXj 

T 7 dy 
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where Y c = (4KT/m e c 2 ) max(r es , r es 2 ) is the Comptonization parameter and 
r es = K es y- Use of equation (7) requires some care since it is meant to describe 
the variation of the radiation temperature when multiple Compton scattering 
is the dominant mechanism to exchange energy between photons and electrons. 
Equation (7) therefore does not apply if either r es < 1 or true emission- 
absorption are important. On the other hand, equation (7) gives the correct 
limit for low optical depth (Y c <C 1 implying T 7 = const), so that one can 
extend the validity of equation (7) to all regimes, provided that no physical 
significance is attached to T 7 where the effective optical depth r e // > 1. 
From equation (1), it follows immediately that 



The boundary conditions for equations (5) and (8) are P = at y = and 
R = R* at y = yi n . We notice that a boundary condition for T 1 must be also 
imposed because the radiation temperature obeys a differential equation; models 
were obtained specifying a value of T 7 at y = 0. The solution to equations (3)- 
(8) provides the variation of L, U, P, T, T 7 and R as functions of the column 
density y and is found numerically; all models refer to a neutron star of R* = 10 
km and M* = 1 M Q . 

We find that two distinct kinds of solutions, "hot" and "cold", always 
exist for any yo provided that the luminosity exceeds a certain limit, which 
depends on yo. The thermal properties of the atmosphere are illustrated in 
figures 1 and 2, where the run of T versus column density is shown for different 
luminosities in the case yo = 20. The "cold" solutions of figure 1 are just those 
already found by ZS and are obtained setting T 7 (y = 0) = [L 00 /(4nR'^a)] 1 / 4 . 
The "hot" solutions of figure 2 exist for values of T 7 satisfying the condition 
T 7 > T crit (Lqo, J/o); here T 7 (y = 0) = 2 x 10 9 K. Temperature is close to T 7 in 
the outer region (y& 23 gem -2 ), while in the dense layers close to y in LTE is 
attained at T ~ 10 7 K. The temperature profile in the hot region resembles that 
one of static atmospheres around X-ray bursting neutron stars (see e.g. London, 



dR 



1 



(8) 



dy 
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Taam & Howard 1986) where the same radiative processes dominate. Cold, 
thermal photons do not propagate outwards because L = for y < y < y in , so 
the hot and the cold zone are thermally decoupled, at least radiatively. On the 
other hand, we have checked that imposing either an ingoing or an outgoing 
flux at yi n does not alter our picture significantly if \L(yi n )\/L 00 ^ 0.1. 

The presence of two possible regimes has a simple interpretation in terms 
of the relative efficiency of the two radiative processes we have considered, 
Compton scattering and free-free emission-absorption. The static atmosphere 
must, in fact, radiate a given luminosity and, for doing that in the scenario 
we are proposing, there are two ways. A first possibility is that a lot of soft 
bremsstrahlung photons are produced in a low-temperature, dense medium 
in which the effective depth is large. This gives rise to a spectrum which is 
essentially blackbody and corresponds to the "cold" solution. Comptonization 
is never dominant because temperature is low and the scattering depth is not 
large enough to make Y c > 1. The "hot" solutions represent the opposite case, 
in which much less energy is generated through bremsstrahlung emission in a 
low-density, hot plasma far from LTE. Comptonization, however, is now so 
efficient that matter and radiation temperatures are everywhere very close and 
the same energy output can be obtained. 

It is possible to get an insight on the existence of high-temperature 
solutions and to give an estimate of the limiting value T cr u by means of 
simple analytical considerations using, for the sake of simplicity, a plane-parallel 
geometry for the atmosphere. For y < y we can safely neglect free-free 
absorption in equation (4) and we get the expressions for L, U and g as functions 
of y and T: 



L = L r 



yo-y 

I 

2/o 



(9) 



W 

u = — 

c 



2y + 3k es y (y - | ) 



(10) 



9 = D V - 



(11) 



where D = GM*m p /2kRl = 8.1 x 10 5 , in c.g.s. units for = 10 6 cm and 
M* = IMq. Neglecting the term KpU, the energy equation becomes a cubic 

equation in x = T 1 / 2 , that can be studied analytically for given values of y 
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and Loo an d treating T 7 as a free parameter. Equation (6) can be cast into the 
form 

x 3 +px + q = (12) 



with 



p=-(T 7 + ^)~-T 7 



B 



U 

y 



U(y) 



where is the total luminosity observed at infinity in units of the Eddington 
luminosity; A = 6.3 x 10 22 and B = 5.1 x 10 25 , again in c.g.s. units. The 
approximated expression for p holds only for T 7 >> 10 8 K, while, in order to 
make the analytical treatment affordable, in the inner part of the atmosphere 
q is set approximately equal to its maximum value, q max = l(yo)- Once yo and 
loo are fixed, equation (12) has one or three real roots, according to the sign of 
the discriminant, and it is easy to prove that all the roots are real only if 

2.5 x 10 8 

(1 + jk es y ) Zoo 

It can also be shown that if just one root is present it is Tj> T 7 , while, when 
condition (13) is satisfied, the three roots have magnitudes T> T 7 , T « T 1 
and T 7 , respectively. The solution T> T 7 is unacceptable since T > T 7 will 
produce a negative radiation temperature gradient (see eq. [7]) and also T « 
T 7 must be discarded because it is inconsistent with our starting assumption 
that absorption could be neglected because the plasma is very hot. Finally, 
the root T^> T 7 is the "hot" solution, which exists only when the radiation 
temperature exceeds the limit given by (13), which represents the analytical 
estimate of T crit . 

In figure 3 we plot the mean energy of the outgoing photons as a function 
of the total luminosity for y = 20. For "hot" solutions only the lower bound 
T crit (dashed line) is shown. As can be seen, while the mean photon energy of 
the "cold" solutions (crosses) monotonically increases with the luminosity, the 
lower bound T crit for "hot" solutions shows the opposite behaviour. Moreover 
the numerical analysis indicates that high temperature solutions may exist only 
for high enough values of /oo and that the critical luminosity , l cr , under which 
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no "hot" solutions exist depends on y : for y = 5, l cr = 2 x 10 2 , while for 
y = 20, l cr = 6 x 10" 3 . 

The anticorrelation between T crit and Z^, and the lack of "hot" solutions 
at low luminosities can be explained as follows. As decreases, W and U 
start to decrease (see eqs. [2] and [10]) and, since Compton heating (which is 
dominant over cooling) becomes progressively smaller in magnitude, the free- 
free emissivity (npaT 4 oc qT 1 ! 2 oc P/T 1 / 2 ) must also decrease for the energy 
equation to be satisfied. This can be achieved only by an increase of temperature 
(and a decrease of density, see eq. [11]) because the pressure profile is nearly 
independent on the thermal properties of the gas. The lack of "hot" solutions at 
small luminosities is due to the fact that, when density is very low, the envelope 
becomes photon starved and even a strong Comptonization is unable to produce 
the required luminosity. 

The existence of both "cold" and "hot" solutions for the same values of 
the flow parameters has been already found in black hole accretion (Park 1990; 
Nobili, Turolla & Zampieri 1991) when both free-free and Compton scattering 
are present. An alternative picture for the production of hard X-rays (~ 100 
keV) from accreting neutron stars has been proposed by Kluzniak & Wilson 
(1991). In their model matter coming from the inner edge of an accretion disk 
hits the stellar surface at a shallow angle, creating a hot equatorial belt in which 
Compton cooling is very efficient. 
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III. DISCUSSION 



The basic limitations of our approach are obviously that the analysis is 
stationary and the spectrum is described just in terms of a mean photon 
energy. The fact that "hot" models may exist only for high enough radiation 
temperatures, T 7 > T cr n ~ 10 9 K at the outer boundary, suggests that, in 
order to get started, an extra energy input, different from that produced by 
the incoming beam, must be supplied, the actual value of T 7 depending on 
the preparation of the system. A physically consistent scenario should be 
investigated in a time dependent picture. Furthermore temperatures are mildly 
relativistic and the expression we used for the Compton heating-cooling term 
is just an approximation. Moreover, pair production-annihilation can not 
be neglected for 10 -2 when the compactness parameter becomes J> 10. 
Above Zqo ~ 0.1, the dynamical effects of radiation pressure and bulk motion 
Comptonization become also important (on this regard see e.g. Zampieri, 
Turolla & Treves 1993 and references therein) and are not included in our 
model. 

Our present results can be summarized as follows. At very low accretion 
rates (the threshold depends on y ) there is a unique solution, the "cold" one 
discovered by ZS (see figure 3). This essentially indicates that the black body 
approximation at very low luminosities, as expected for instance in isolated 
neutron stars accreting the interstellar medium is reasonable (see e.g. Ostriker, 
Rees & Silk 1970; Treves & Colpi 1991; Blaes & Madau 1993). The "hot" 
and "cold" solutions coexist for a certain range of luminosities and the mean 
photon energies of the two modes approach each other for increasing l^. The 
"cold" solution becomes hotter for increasing accretion rate while the "hot" one 
softens and this behaviour is opposite to the one exhibited by the hot, shocked 
solutions of Shapiro & Salpeter. 

Even in the absence of frequency-dependent calculations it is tempting to 
associate the "cold" solution with the spectral states observed in the Soft X- 
Ray Transient source Aql X-l, either in outburst or in quiescence (Verbunt 
et al. 1993). The existence of "hot" solutions for intermediate luminosities 
with temperatures ~ 100 keV suggests that a class of hard X-ray sources, 
which possibly have not yet been discovered, may exist. The transition between 
the "hot" and the "cold" regime, even at luminosities where the two solutions 
are rather different, may be expected in a time-dependent scenario. Non- 
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stationary calculations are also needed to explore the stability properties of the 
two solutions. 

In conclusion the discovery of "hot" solutions in spherical accretion onto 
neutron stars deserves further theoretical consideration, a program on which we 
are presently actively working. 
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FIGURE CAPTIONS 

Figure 1. Temperature vs. column density of "cold" solutions for = 7x 10~ 3 
(continous line), 1^ = 2 x 10~ 2 (dashed line) and = 7 x 10~ 2 (dashed- 
dotted line); here y = 20. 

Figure 2. Same as in figure 1 for "hot" solutions. 

Figure 3. Mean energy of the outgoing photons vs. total emitted luminosity 
for y = 20; crosses refer to "cold" models. The dashed line represents the 
lower limit for the existence of "hot" solutions given by equation (13). 
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